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A  CON8ERVATION-OF- VELOCITY  LAW  FOR  INVISCID  FLUIDS 


1.  Introduction 

Conservation  laws  often  provide  keys  to  the  understanding  of  basic 
phenomena  throughout  physics.  In  the  study  of  water  waves  a  single 
conservation  law  -  energy  -  is  sufficient  to  derive,  for  example.  Green's  lav 
for  long  linear  waves  (Green,  1838)  and  an  analogous  law  for  solitary  waves 
in  a  channel  of  gradually  varying  section  (see  Miles,  1979)*  Lamb  (1932) 
invokes  two  linear  conservation  laws  (mass  and  velocity  or  momentum)  in  an 
approximate  way  to  derive  reflection  coefficients  for  long  linear  waves 
encountering  a  barrier.  The  results  agree  well  with  the  more  general  linear 
solution  of  Bartholomeusz  (1958)  when  the  waves  are  long  enough.  Witting 
(1981)  invokes  the  two  linear  lavs  and  a  quadratic  law  to  derive  reflection 
and  re-reflection  coefficients  for  long  linear  waves.  He  shows  that  failure 
to  invoke  the  velocity  conservation  lav  (while  retaining  mass  conservation) 
removes  the  re-reflected  wave,  and  alters  the  properties  of  the  reflected 
wave. 

In  this  brief  report  we  show  that  the  "conservation  of  velocity"  law  as 
used  by  Witting  is  a  manifestation  of  a  very  general  linear  conservation  law 
for  invlscid  water  waves.  For  irrotational  flows  it  is  a  variant  of 
Bernoulli's  lav,  but  the  general  form  holds  for  rotational  flows  as  well.  We 
believe  that  this  conservation  of  velocity  lav  nay  be  novel  (at  least  in  its 
present  form).  In  Section  2  we  derive  the  lav  and  discuss  its  connection  to 
Bernoulli's  law.  Some  applications  are  discussed  in  Section  3* 


2.  Equation*  of  Motion 
A.  General  Derivation 

We  choose  a  coordinate  system  as  illustrated  in  Fig.  1.  with  x 
horizontal,  y  vertical,  and  n  (x,t)  the  elevation  of  the  surface  (a  free 
surface  or  an  arbitrary  continuous  surface  of  markers  moving  with  the  fluid) 
above  a  line  y constant.  The  Euler  equations  for  the  fluid  and  surface 
elevation  are 
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We  now  define  x  and  t  derivatives  at  a  fixed  value  of  x,  but  moving 
vertically  vith  the  surface.  For  arbitrary  f(x,y,t),  let  the  derivatives 
d/dx  and  d/dt  be  defined  as  follows: 
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From  this  point  on,  derivatives  d/dx  and  d/dt  will  illicitly  take  y»n,  as  in 
(4)  and  (5). 
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At  the  surface,  (l)  -  (3)  become 
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Invoking 


3  3n  _  3  3n  ,  (0\ 
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and  p  *  constant,  (6)  -  (8)  give  straightforwardly 

ft <u  +  4ir v) "  -  (S^  *ur£*f*  «")•  <10> 

In  the  more  general  case  p  *  p(p),  p/P  in  (10)  becomes  /dp/p.  The  general 
result  depends  neither  on  irrotationality  nor  inconpressibility.  The  operand 
of  the  time  derivative  in  (10)  can  be  expressed  as  follows.  I«t  Vip  denote 
the  tangential  coqponent  of  the  surface  velocity. 

Then 

u  +  v  ■  vT  sec  8,  (11) 

with  tan  8  the  surface  slope  (see  Figure  1).  Thus  the  left  hand  side  of  (ll) 
is  proportional  to  the  tangential  surface  velocity. 

B.  Connection  to  Bernoulli’s  Lav 

For  irrotational  flows  Bernoulli's  Lav  is: 


=  -  l§  (u2  +  v2)  +  *  +  gy  +  F(t)| 


(3  2) 

where  ♦(x,y,t)  is  the  velocity  potential  and  F(t)  is  an  arbitrary  function  of 
the  time.  By  using  (5)»  (12)  becomes 

5  IT  -  -  l|  (“2  ♦  \  ♦  ST  ♦  «t>!  ♦  V  |2.  (13) 

Now  take  d/dx  of  (13)  and  interchange  d/dx  and  d/dt: 


ft  (dVdx)  -  *  fj  ^2  +  v2)  +  £  ♦  gn  -  v  j£).  (1U) 

The  kinematic  condition  along  y  ■  n  is: 

v(x,n,t)  =  v  ■  Bn/at  +  u  Jn/ix,  (15) 
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so  that  the  substitution  of  (15)  for  9n/at  in  (lU)  gives: 
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Using  (U),  we  interpret  d+8/dx  as  u8  ♦  (3n/3x)va,  and  so  (16)  is  Just 
(10).  Consequently,  we  may  Interpret  (10)  as  a  Bernoulli  Law.  Note, 
however,  that  the  starting  point  (12)  is  valid  only  for  irrotat tonal  flows, 
while  the  result  (10)  is  valid  for  rotational  flow  as  well.  Of  course,  the 
starting  point  for  (10)  is  the  set  of  Euler  Equations,  and  so  (10)  is 
strictly  valid  only  for  flows  in  inviscid  fluids. 


3*  Applications  and  Discussion 


A*  Sample  Calculations  using  Conservation  of  Velocity. 

The  velocity  conservation  lav  represented  by  (10)  is  one  of  two  general 
linear  conservaton  lavs  for  tvo  dimensional  lnviscid  flows.  The  other  is  the 


well  known  conservation  of  mass: 

St  *  5T  ■  0  <1T> 

where  A  is  the  cross-sectional  area  above  still  water  level  /ndz,  where  z  is 
the  coordinate  perpendicular  to  the  plane  of  fluid  motions,  and  F  is  the 
volume  flux  //  u  dy  dz,  taken  over  the  whole  fluid  at  a  fixed  value  of  x. 
Various  forms  of  the  Bousslnesq  Equations  for  irrotational  water  waves 
conserve  two  linear  quantities,  corresponding  to  the  mass  of  (IT)  and  the 
velocity  of  (10). 

Simplifications  to  the  Bousslnesq  Equations  that  permit  waves  to  travel 
in  one  direction  only,  such  as  that  derived  by  Kbrteweg  and  de  Vries  (1895), 
conserve  neither  the  mass  of  (17)  nor  the  velocity  of  (10).  Bousslnesq 
(1872),  Miles  (1979)*  and  a  few  researchers  in  the  intervening  centuxy 
recognized  that  mass  was  not  conserved  in  monodirectlonal  formulations,  and 
added  a  reflected  wave  to  balance  a  mass  budget.  The  energetic  part  of  the 
disturbance  in  the  raonodirectional  formulation  (a  modified  KdeV  Equation  in 
Miles's  work)  was  left  alone. 

Witting  (1982)  has  developed  a  model  that  incorporates  mass  and  velocity 
conservation  explicitly.  The  model  equations  of  motion  are: 
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The  system  of  equations  (l8)-(2l)  is  closed  by  relating  u  to  q  through  power 
series  expansion.  Equation  (18)  is  a  restatement  of  (17),  and  (19)  follows 
from  (10)  and  (20).  Some  calculations  using  the  theory  and  numerica1  methods 
described  by  Witting  (1982)  were  designed  to  model  experiments  of  Chang,  et 
al.  (1979)  in  a  converging/diverging  channel. 

Figure  2  displays  a  profile  n  ■  n  (x,tQ)  of  a  wave  in  a  converging 
j  channel.  In  the  figure  we  mark  the  crest  location  by  x  ■  0,  and  express  n 

and  x  in  units  of  still  water  depth  hQ.  The  particular  time  tQ  is  chosen 
to  match  the  time  at  which  our  crest  passed  the  location  of  a  wave  probe  of 
Chang,  et  al. ,  who  measure  n  *  n  (xQ,t).  The  probe  record  is  chosen  from 
Chang,  et  al.  's  Figure  7c.  The  profile  of  n  in  Figure  2  resembles  the 
experimental  record.  Both  are  asymmetrical,  about  to  the  same  degree,  if 
allowance  is  made  for  Chang's  slightly  noisy  records  and  for  the  small 
differences  that  must  arise  in  comparing  a  space  record  with  a  time  record. 

In  order  to  identify  the  source  of  the  asymmetry,  we  also  plot  q 
(nondimens ionall zed  by  (ghQ)1/2l,  (n+q)/2,  and  5(n-q)  in  Fig.  2.  If  all 
disturbances  were  linear  long  waves,  (n+q)/2  would  represent  a  disturbance 
traveling  to  the  right  (n  in  phase  with  q),  and  5(n-q)  would  represent  a 
disturbance  traveling  to  the  left  (n  in  phase  with  -q,  and  scale  expanded 
tenfold).  In  Fig.  2  the  (n+q)/2  profile  is  still  slightly  asymmetrical, 
which  agrees  with  the  theoretical  profile  that  Chang,  et  al.  derive  from  a 
Korteweg-deVries  theory  in  comparing  their  experimental  data  to  theory.  The 
5(  n-q)  profile  clearly  shows  a  reflected  wave,  which  the  Korteweg-deVries 
1  theory  does  not  describe.  (We  should  note  that  only  for 

x<-U  is  the  ( h+q) /2  profile  small  enough  that  5( n-q)  unambiguously  signals  a 
reflected  wave.  For  x>-4  a  decomposition  into  right-  and  left -going  waves 
needs  a  more  detailed  theory  than  linear  long  wave  theory.) 

If  mass  conservation  is  added  as  in  Miles  (1979) »  the  reflected  wave  will  be 
taken  into  account,  perhaps  adequately. 
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DIMENSIONLESS  DISTANCE  FROM  CREST 

Figure  2  -  Solitary  wave  in  a  converging  channel.  The  geometry  is  set  to 
that  of  of  Chang,  et  al.  (1979).  which  has  straight  vails  converging  from  the 
end  of  an  entrance  section  at  an  angle  of  1.07**°.  The  elevation  above 
still  water  level  and  distance  from  the  crest  are  measured  in  units  of  hQ, 
the  initial  depth  (■O.UOm).  The  velocity  variable  q=  ug+n'vg  is  measured 
in  units  of  (ghQ)l/2.  At  the  crest  the  channel  is  O.U3*»  as  broad  as  the 
entrance  setion,  where  the  channel  walls  are  parallel.  Initial  conditions 
were  set  to  produce  a  solitary  wave  in  the  entrance  section  having 
nondimensional  amplitude  0.197.  The  crest  in  the  figure  has  nondimensional 


elevation  0.321 


Figure  3  displays  profiles  of  n,  q,  (h+q)/2  and  5(  n-q)  in  a  diverging 
channel*  Again  x  *  0  narks  the  crest  location.  The  time  at  which  the 
profile  "snapshots"  are  taken  matches  the  time  that  the  crest  passed  one  of 
the  wave  probes  of  Chang,  et  al.  (the  probe  record  is  shown  in  the  fifth 
display  of  their  Figure  6).  The  computed  profile  of  n  in  Figure  3  again 
reproduces  the  essential  features  of  the  experimental  record.  The  wave  is 
asymmetrical  in  the  oppsite  sense  of  Figure  2,  and  has  an  oscillatory  tail 
that  is  mostly  negative.  By  contrast,  the  Korteweg-deVries  theory  of  Chang, 
et  al. ,  which  conserves  neither  mass  nor  velocity,  gives  an  almost 
symmetrical  solitary  wave  with  no  reflected  wave  and  no  tail.  From  the 
bottom  curve  of  Figure  3  we  see  that  some  of  the  disturbance  behind  the 
solitary  wave  is  a  reflected  wave  of  depression.  From  the  ( n+q) /2  curve  we 
see  that  the  rest  of  the  disturbance  behind  the  solitary  wave  is  an 
oscillatory  wave,  mostly  negative,  that  is  traveling  toward  the  right,  i.e. 
along  with  the  solitary  wave. 

Thus,  the  calculations  displayed  in  Fig.  3,  which  conserve  both  mass  and 
velocity,  reproduce  the  features  behind  the  solitary  wave  that  show  up  in  the 
laboratory  experiments.  These  features  appear  to  be  the  result  of  waves 
traveling  both  to  the  left  and  to  the  right.  We  believe  that  the  rightward 
traveling  wave  is  composed  of  both  a  modified  solitary  wave  and  a  re¬ 
reflected  tail  akin  to  that  described  by  Witting  (19BI). 

B.  Discussion 

The  principal  new  result  of  this  research  is  the  conservation  of  velocity 
law  derived  under  general  circumstances  (Eq.  10).  It  can  be  used  explicitly 
in  numerical  calculations,  or  can  be  used  as  a  test  for  the  accuracy  of 
calculations  that  do  not  make  use  of  it  explicitly.  Because  it  is  embedded 
in  the  Euler  Equations,  direct  use  of  these  equations  should  conserve  it 
automatically.  Certain  approximations,  such  as  the  modified  Korteweg-deVries 
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DIMENSIONLESS  DISTANCE  FROM  CREST 

Figure  3  -  Solitary  wave  in  a  diverging  channel.  The  geometry  is  set  to  that 
of  Chang,  et  al. ,  (1979),  which  has  straight  walls  diverging  from  the  end  of 
an  entrance  section  at  an  angle  of  1.07^°.  As  in  Fig.  2,  hQ  and 
(ghQ)1/2  are  the  units  of  length  and  velocity  (but  with  hQ  ■  0.30m  here). 

At  the  crest  the  channel  is  6.822  times  broader  than  the  entrance  section. 

In  the  entrance  section  the  distrubance  was  a  solitary  wave  of  amplitude 
0.386.  The  crest  in  the  figure  has  non dimensional  elevation  0.109* 


Equation  for  wave  propagation  in  a  channel  of  varying  section,  violate 
conservation  of  mass  and  velocity.  A  model  intended  to  account  for  a  broader 
range  of  phenomena  as  described  above  should  conserve  both  mass  and  velocity. 
Equation  (10)  with  (IT )  provides  an  explicit  (though  not  unique)  mechanism 
for  performing  this  task. 
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